Mon. Not. R. Astron. Soc. 000, dini (2003) Printed 2 February 2008 (MN WF^ style file v2.2) 



Evidence for a multi— zone warm absorber in the 
XMM— Newton spectrum of Markarian 304 



E. Piconcelli^*, E. Jimenez-Bailon^, M. Guainazzi^, N. Schartel^, 
P.M. Rodriguez-Pascual^, M. Santos-Lleo^ 

1 XMM-Newton Science Operation Center/RSSD-ESA, Apartado 50727, Madrid 28080, Spain 



Accepted. Received 



ABSTRACT 

We present a XMM-Newton observation of Markarian 304, a Seyfert 1 galaxy at z = 
0.066. The EPIC data show that MKN 304 is affected by heavy (Nr « lO^^ cm-^) 
obscuration due to ionized gas. A two-phase warm absorber provides an adequate 
parameterization of this gas. The ionization parameter of the two components is ^ « 6 
erg cm~^ s~^ and ^ « 90 erg cm~^ s~^, respectively. The observed continuum photon 
index (F « 1.9) is typical for Seyfert 1 galaxies. Two significant emission lines are 
detected at 0.57 keV and 6.4 keV, respectively. The former is mostly likely due to He- 
like oxygen triplet emission arising from an ionized plasma (maybe the warm absorber 
itself) . The latter is due to fluorescent emission of K-shell iron in a low-ionization state 
(Fel-XV). The upper limit for the line width of aKa < 0.18 keV most likely rules out 
an origin in the inner parts of the accretion disk. Interestingly, the strength of such line 
is consistent with the possibility that the emission is produced in the warm absorber 
itself. However, a substantial contribution from the torus is plausible too. We have 
also found a weak (4% of the primary continuum) soft excess emission component. 
The presence of this excess could be explained by either emission/scattering from a 
warm gas or partial covering, or a combination of them. 

Key words: galaxies: active - galaxies: individual: Markarian 304 - galaxies: Seyfert 
- X-rays: galaxies. 



1 INTRODUCTION 

The bulk of the X-ray emission in active galactic nuclei 
(AGNs) is produced via Comptonization of UV photons in 
the inner parts of the accretion flow close to the central su- 
permassive black hole. The resulting power law spectrum 
is modified by a number absorption and emission spectral 
features due to the reprocessing of this primary continuum 
(Mushotzky, Done & Pounds 1993). They are clear signa- 
tures of the presence of large amounts of gas characterized 
by different physical properties in the circumnuclear region. 

One of the most common spectral features is the so- 
called warm (i.e. partially ionized) absorber (Halpern 1984; 
Pan et al. 1990, Turner et al. 1993). Such a component has 
been found to contribute significantly to the opacity of the 
X-ray primary emission in ~ 50% of the Seyfert 1 galaxies 
with column densities up to 10'^'^ cm~^ (Reynolds 1997). 
In particular, OVII (0.739 keV) and OVIII (0.871 keV) ab- 
sorption edges resulted to be the most prominent absorption 
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signatures of the ionized gas in low-resolution ASCA obser- 
vations. 

The advent of grating spectrometers on-board XMM- 
Newton and Ctiandra has dramatically enriched our knowl- 
edge in this field. High-resolution soft X-ray spectra of 
bright Type 1 AGNs present a wealth of absorption lines 
(e.g. Collinge et al. 2001; Kaspi et al. 2001; Kaastra et al. 
2002), which allow to infer a few remarkable conclusions: (i) 
warm absorbers usually show a very complex structure, with 
multiple regions characterized by different ionization states, 
column densities and velocities; (ii) the strongest features 
are due to He- and H-like states of O, Ne, Mg and Si; in 
addition a blend of M-shell iron inner shell transition (un- 
resolved transition array or UTA) has been often detected 
in the range « 0.70-0.80 keV; {iii) the observed blueshift of 
the absorption lines imply they are originated in an outflow 
with mean velocities in the range « 100-2500 km/s. 

However, two important issues regarding the warm ab- 
sorbers still remain open, i.e.: the exact location of the ab- 
sorbing medium and the physical connection between the 
absorbing materials in the UV and X-ray band. A likely 
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range of distances of the ionized absorbing medium from 
the central source is r w 0.01-10 pc (KroUk 2002; Behar et 
al. 2003). 

By correlating HST and ASCA data Crenshaw et al. 
(1999) found that all Scyfcrt galaxies with evidence of UV 
absorption also exhibit X-ray ionized absorption features. 
Furthermore, some sources in which UV and X-ray absorb- 
ing components share the same outflow velocity space have 
been observed (Blustin et al. 2003; Collinge et al. 2001; 
Kaastra et al. 2002). However, no cloarcut identification of 
an X-ray warm absorber with an UV counterpart has been 
reported yet (see Crenshaw, Kraemer & George 2003 for a 
review) . 

Unfortunately, the sensitivity of current grating spec- 
trometers allows us to analyze just a handful of very bright 
X-ray sources. Therefore CCD cameras spectroscopy still 
represents the best available tool to deduce information 

about the physical properties of the gas in the immediate 
surroundings of the central X-ray emitting region for most 
AGNs (i.e. those with a F2-10 ^ 1 mCrab^). We present 
here the analysis of an EPIC observation of Markarian 304 
(F2-10 « 0.1 mCrab). 

Markarian 304 (MKN 304, PG2214+139) is a Scyfcrt 
galaxy at 2 = 0.066. Due to the presence of broad emission 
lines in its optical spectrum, MKN 304 has been classified 
as a Type 1 object (Osterbrock 1977). lUE spectroscopy of 
this source revealed the presence of strong bump in the UV 
(Clavel & Joly 1984), without absorption lines. On the con- 
trary, a recent FUSE observation of MKN 304 (Kriss 2002) 
disclosed broad absorbing features (with a velocity spread 
of 1500 km/s) due to multiple OVI absorption components 
blended together. 

MKN 304 was discovered as X-ray emitting source by 
Uhuru (Tananbaum et al. 1978). Later X-ray observations 
of this source were performed by Einstein (Kriss, Canizares 
& Ricker 1980) and ROSAT (Rachen et al. 1996): both re- 
ported evidences of heavy obscuration, with a very flat slope 
(r - 0.07) in the soft X-ray band. Surprisingly, MKN 304 
was not targeted by ASCA and BeppoSAX and, therefore, no 
further study of its unusual flat continuum was performed. 

Throughout this paper we assume a flat A CDM cosmol- 
ogy with (Omj^a) = (0.3,0.7) and a Hubble constant of 70 
km s"^ Mpc"^ (Bennett et al. 2003). 



2 XMM-Newton OBSERVATIONS AND DATA 
REDUCTION 

MKN 304 was observed by the European Photon Imaging 
Camera {EPIC) on board the XMM-Newton space obser- 
vatory (Jansen et al. 2001) between 2002 May 12 15:46:29 
(UT) and 2002 May 13 00:56:34 (UT). The observation was 
performed with the EPIC PN camera (Struder et al. 2001) 
operating in Full Window Extended mode and with the 
EPIC MOS cameras (Turner et al. 2001) operating in Full 
Window mode. All the three EPIC cameras were equipped 
with the "Thin" blocking optical filters. 

^ 1 Crab corresponds to a 2-10 keV flux of ~ 2 x 10~* erg cm~^ 
s-i. 



The data files were reduced and analyzed using the 
standard Science Analysis System {SAS) v.5.4.1. We used 
the EPCHAIN and EMCHAIN tasks for processing of the 
raw PN and MOS data files, respectively, in order to gen- 
erate the relative linearized event files. X-ray events cor- 
responding to patterns 0-12 (i.e. single-, double-, triple- 
, and quadruple-pixel events) for the MOS cameras were 
selected; for the PN, patterns 0-4 events (i.e. single- and 
double-pixel events) were used. The subsequent event selec- 
tion was carried out taking into account the most updated 
cahbration files at the time the reduction was performed 
(September 2003) . All known flickering and bad pixels were 
removed. The event lists were furthermore filtered to ignore 
periods of high background fiaring (see Appendix A). Af- 
ter this screening process (and taking into account satellite 
dead times) the efi'ective net exposure times for the MOSl, 
M0S2 and PN detectors are 33.4, 35.5 and 29 ks, respec- 
tively. No exposure is affected by pile-up. The peak of the 
X-ray emission is located at 02000 ~ 22''17'"12.1° and S-zom 
— -1-14° 14'21.7", i.e. at a distance of 1.5 arcsccs from the 
coordinates of the optical nucleus of MKN 304 (Clements 
1981). This is consistent with typical attitude reconstruc- 
tion uncertainties. 

Since the cleaned RGS data did not yield enough pho- 
tons for a meaningful spectral analysis, they are not dis- 
cussed in the following. 

The source spectra were extracted from the final fil- 
tered event lists using a circular region of 24 and 21.5 arc- 
sees radius for PN and MOS, respectively, centered on the 
observed position of MKN 304. Backgrounds were estimated 
from source-free similar regions close to the source on the 
same CCD. Response matrices and ancillary response files 
were generated using the RMFGEN and ARFGEN tools in 
the SAS software. Since difference between the MOSl and 
MOS 2 response matrices are a few percent, we created a 
combined MOS spectrum and response matrix. The PN and 
MOS spectra were then fitted simultaneously. Both spectra 
were grouped to give a unnimum of 40 counts per bin in or- 
der to apply statistics. Given the current calibration un- 
certainties and the detector sensitivities, events outside the 
0.3-12 keV range were ignored in the PN spectrum, while, 
for the MOS, we retained the 0.6-10 keV band. 



3 X-RAY SPECTRAL ANALYSIS 

Spectral analysis weis performed using XSPEC v. 11.2 (Ar- 
naud 1996). All the models presented in this Section and 
in the following one include absorption due to the line-of- 
sight Galactic column of Nh — 4.96 x 10^" cm^^ (Dickey 
& Lockman 1990). The photoelectric absorption cross sec- 
tions were taken from Morrison & McCammon (1983). The 
quoted errors refer to the 90% confidence level for one inter- 
esting parameter (i.e. Ax^ = 2.71; Avni 1976). 

The 0.5-2 keV and 2-10 keV PN lightcurves were vi- 
sually inspected for checking variability. We did not detect 
any significant change in both count rates during the w 9 hr 
observation. 

We first fitted the EPIC spectra with a power law model 
in the 2-12 keV band. This simple model yielded a very poor 
fit {x^{i')= 469(256)) and an extremely flat photon index of 
F « 0.95. Indeed the convex shape of the data-to-model 
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Figure 1. Upper panel: Continuum power law fit to the 2—12 IceV Figure 2. Upper panel: The PN (top) and MOS (bottom) spec- 

band of the PN (top) and MOS (bottom) spectra of MKN 304, ex- tra of MKN 304 fitted by a single warm absorber model (model 

trapolated over 0.3—12 (0.6—12) keV for PN{MOS). Bottom panel: B in Table 1). Bottom panel: Residuals in units of standard devi- 

Residuals in units of standard deviations. ations. They clearly indicate that this model cannot account for 

the complex spectrum of MKN 304. 



residuals suggests the presence of heavy obscuration. Fig.Q 
displays the results of this fit extrapolated over the 0.3- 
12 keV bandpass. The large deficit in the 0.5-2 keV band 
clearly indicates the presence of a strong absorption, while 
the systematic positive residuals at ^ 0.6 keV may reveal 
an emission component in excess of the underlying absorbed 
continuum. Fitting the data over the 0.3-12 keV band with 
a power law yielded a X^iv) ~ 2483(367) (see model A in 
Table EJ. 

We thus accounted for this ionized absorber component 
with the multiplicative absori model in XSPEC (Done et 
al. 1992). In this model the state of the warm absorber is a 
function of the ionization parameter ^ defined as ^ = L/nr^, 
where L is the isotropic luminosity of the ionizing source in 
the interval 5 eV to 300 keV, n is the number density of 
the warm plasma and r is the distance between this latter 
and the central source. In the fit solar elemental abundances 
were assumed and ^ was left as a free parameter. On the con- 
trary, we fixed the temperature of the absorber to T = 1.5 
X lO'^ K since the statistics did not allow to simultaneously 
and effectively constrain T and ^. Such a value of T falls in 
the typical range observed for a warm absorber (Reynolds 
& Fabian 1995; Krolik 2002 and reference therein). We also 
performed a further set of fits assuming values of tempera- 
ture in the range 5 X 10'' < T < 10^ K, from which turned 
out that the ionization parameter remains almost constant 
(i.e. ^ « 90-100 erg cm~^ s"'^) to the variations of temper- 
ature. This fit (indicated as model B in Table resulted in 
a dramatic statistical improvement i.e. at > 99.99% confi- 
dence level according to an F-test once compared to model 
A, with a X^{v) ~ 617(365). Nevertheless, the value of /v 
~ 1.7 remains unacceptable and a large absorption feature is 
evident in the soft portion of the spectrum (see Fig. 2). The 
breadth and depth of this feature suggest the existence of 
another ionized absorption component with a different value 
of ^. Furthermore the derived slope F = 1.40±0.04 is signif- 
icantly fiatter than usually observed in unabsorbed Seyfert 
galaxies (Malizia et al. 2003) and quasars (Piconcelli et al. 
2003), i.e. (F) ^ 1.9. 

We then added another absori component in the spec- 
tral model. First we fitted the spectrum fixing the tempera- 
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Figure 3. Left panel: contours at the 68%, 90% and 99% confi- 
dence level for energy and flux of the He-like oxygen line. Right 
panel: contours at the 68%, 90% and 99% confidence level for 
energy and flux of the iron line. Both fluxes are in units of 
photons / cm^ /s. 



ture of this new component to the same value of the previous 
one, i.e. T — 1.5 x 10^ K. The fit gave a minimum X^{v) = 
477(363), with a statistical improvement significant at > 
99.99 confidence level. However, an even better value [x^ 
— 441) was obtained fixing the temperature of the second 
ionized component to a lower value, in particular with T 
= 3 X 10'* K. The resulting spectral parameters of this fit 
(Model C) are also listed in Tabled The value was re- 
duced by Ax^ = 176 respect to Model B, thus the addition 
of the low-temperature ionized absorption term is significant 
at a level larger then 99.99% confidence level (F-test value 
= 72). 

However the reported value of X^/'^ = 1-22 (null hy- 
pothesis probability 0.003) is not yet very satisfactory and 
positive systematic residuals are still present in the data- 
to-model ratio below 1 keV. 

We added a second power law accounting for this addi- 
tional soft X-ray emitting component. The relative photon 
index was constrained to be that of the primary continuum 
but the normalization was left free to vary. 

We obtained again a very significant improvement in 
the fit quality at the > 99.99% confidence level {F-test = 
66) with reduced x^ = 373 for 362 degrees of freedom (see 
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Table 1. X-ray spectral model fitting results. Parameters marked with T5(T4) are referred to the high-(low-) temperature component of 
the warm absorber. In model E the photon index of the unabsorbcd power law was fixed to the value of the primary absorbed continuum. 
See text for details. 



model D in Table 1). The resulting photon index is F = 
1.92±0.08 and the normalization of the unabsorbed power 
law is « 4% of the absorbed power law component. The best 
fit column densities of the two ionized absorbers increased 
to Ng* = 1.3±0.3 X 1022 cm-2 j^ts ^ io.9«;^ x IO22 
cm-2 £|_^j. ^j^g lukewarm (i.e. with T = 3 x 10* K) and hot 
(i.e. with T = 1.5 x 10^ K) component, respectively. The 
ionization parameters turn out to be ^"^^ — 3.8+2-4 erg cm-2 
s"^ and = llO.St? * erg cm-2 ^-i 

We attempted an alternative fit modelling the soft ex- 
cess with a blackbody component attenuated by both ion- 
ized absorbers, and with all fit parameters free to vary, which 
yielded a minimum = 399(361) and, therefore, signif- 

icantly worse than model D. We also tried to substitute the 
"lukewarm" ionized absorbing component in model D with 
a neutral one, but the quality of the fit resulted very poor 
with a final X^iy) ~ 532(363) (null hypothesis probability 
1.5 X 10"®). 

Model D provides an excellent and physically plausible 
parametrization of the overall EPIC spectrum of MKN 304, 
with a null hypothesis probability equal to 0.5. Visual in- 
spection of the data-to-model ratio residuals suggests the 
presence of two emission features, respectively, at ~ 6 keV 
and ~ 0.5 keV (observer-frame). The former is likely due 
to the fiuorescence emission from the iron K-shell as com- 
monly observed in most of Seyfert galaxies (Nandra et al. 
1997). The latter is likely due to the helium-like oxygen 
triplet of emission lines centered at 0.57 keV (rest-frame) 
as found in other Seyfert galaxies (e.g. Kaastra et al. 2002) . 
Thus the data led us to add two narrow gaussian lines in 
model D to account for them. The best fit rest-frame energy 
of the Fe Ka hne is 6.34±0.06 keV with an equivalent width 
of EW = 79±35 eV. This hne is significant at 99.9% confi- 
dence level (Ax^ = 13 for two additional parameters; _F-test 
value = 7.3), the centroid energy corresponds to low ioniza- 
tion states, i.e. Fel-XV (Makishima 1986). We obtained an 
upper limit for the line width of ukc < 0.18 keV. The rest- 
frame position of the soft X-ray line {E = 0.57lg Qi keV) is 
fully consistent with OVII emission. Given the reduction of 
of Ax^ = 25 for two additional degrees of freedom, the 
detection of this line is significant at > 99.99% confidence 
level based on an _F-test. In Fig. |5] are shown the 68.3%, 
90% and 99% confidence levels for two interesting parame- 
ters between the line energy and the fiux of both emission 
features assuming model E in Table 1. 

Fig-IUshows the PN and MOS spectra and the spectral 
fit which includes the two gaussian emission line. We con- 
sider this spectral model (Model E in Table l)as the best fit 




channel energy (keV) 

Figure 4. Unfolded best-fitting _BP/C spectrum of MKN 304 fit- 
ted with two power laws plus a double warm absorber component 
and two gaussian lines (model E in Table 1). The photon index 
of the unabsorbed power law is fixed to the value of the primary 
absorbed continuum (see text for details). 



to the EPIC data of MKN 304. This fit yields a minimum 
X^{v) = 335(358) with a null hypothesis probability of 0.82. 

Once Model E was assumed, we measured a fiux of 
Fo.5-2 = 3.06 X 10"" erg cm"2 s'^ and ^2-10 = 3.26 
X 10"^2 ej-g g-i ^jjg (0.5-2 keV) and high 

(2-10 keV) energy bands, respectively. After correction for 
both Galactic absorption and warm absorber column den- 
sities, these correspond to luminosities of Lo.5-2 = 3.9 x 
10*'^ erg s"^ and L2-10 = 4.8 x 10*^ erg s"^, respectively. 
During the XMM-Newton observation the soft X-ray fiux 
of MKN 304 was therefore a factor of 0.5 lower than mea- 
sured2 by Einstein on 1979 and a factor of 1.4 higher than 
the value obtained by ROSAT on 1992. Unfortunately, due 
to the lack of spectral information we cannot investigate the 
possibility that these variations are linked to changes in the 
column density of the observer or, alternatively, caused by 
the activity cycle of the X-ray continuum source. 



2 In order to compare the fluxes measured by different satellites 
we assumed a power law model with F = 1.8 and Nh= N^'''. 
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4 DISCUSSION 

The spectrum of MKN 304 unveiled by the XMM- 
Newton observation discussed in this paper is complex both 
in the soft and in the hard energy band (Fig. 1). 

Our analysis fully confirms previous results on nar- 
rower energy ranges which reported a very flat continuum 
for this source, and suggested the existence of heavy obscu- 
ration. However, due to the limited bandpass and the poor 
sensitivity of the previous observations, it was impossible 
to accurately investigate the properties of the absorption, 
and discriminate between a cold (neutral) or warm (ion- 
ized) medium. This XMM-Newton observation of MKN 304 
has clearly revealed that the strong obscuration (Nh ~ 10^^ 
cm~^) occurs in an ionized gas. 

We have indeed found that the best fit model to the 
EPIC data consists of a two-phase warm absorber with a 
high- and a low-ionized component, a faint soft excess com- 
ponent and two emission lines at ~ 0.57 keV and ~ 6.4 keV, 
respectively (see Model E in Table 1). 



4.1 Warm absorber 

We fixed a priori the temperature values at T = 1.5 x 10''' 
K and T = 3.0 x lO" K for the two "phases" of the ion- 
ized plasma since they represent typical values observed in 
sources with a multi-phase warm absorber (Netzer et al. 
2003; Kaastra et al. 2002). From the spectral analysis we 
derived i^'' = 89.311^:^5 erg cm'^ s"' and = 5.9+^ * erg 
cm~^ s~^ for the hot and the cold component, respectively 
(see model E in Table 1). 

Our findings for MKN 304 appear to be in agreement 
with a "multi-zone" structure for the ionized absorber in 
Seyfert galaxies as discovered by high resolution spectro- 
scopic observations recently performed with grating spec- 
trometers on-board of XMM-Newton and Chandra. They 
indeed resolved the complex structure of the AGN circum- 
nuclear absorbing material into a multi-phase plasma with 
several orders of magnitude spread in ^ and systematic 
blueshifts (Blustin et al. 2003; Kaastra et al. 2002). This 
observational evidence leads to interpret the warm absorber 
phenomenon essentially in terms of a photoionized outflow- 
ing wind consisting of different inhomogeneous co-existing 
regions (or shells) characterized by a broad distribution in 
5, temperature and density (Krolik 2002; Behar et al. 2003; 
Elvis 2000). 

Interestingly, while our best-fit values for ionization 
parameters of the two absorbing components are similar 
to those found in other Seyfert 1 galaxies (Kaastra et al. 
2002; Blustin et al. 2003), the corresponding column densi- 
ties (N5^ « 9 X 10^^ cm-2 and Ng* fa 2 x 10^^ cm'^) are 
amongst the highest seen by Chandra and XMM-Newton so 
far. 

Given these high values of Nh, photoionization edges 
are expected to dominate the absorption features in 
MKN 304. In fact, at high column densities the UTA ab- 
sorption lines saturate. Edges (alongside with 2p — nd (with 
n > 3) absorption lines) provide the major contribution 
to the absorption spectrum (Behar, Sako & Kahn 2001). 
Furthermore as shown in Fig. the absorption feature in 



MKN 304 is remarkably broad and deep and, hence, hard to 
reconcile with what simply expected from UTA^. 

Besides absorption, also emission lines are expected to 
arise from the warm medium. Due to the underlying illumi- 
nating bright continuum and the forest of absorption lines 
such emission features are however difficult to detect in low- 
resolution data. In particular He-like triplets of O and Ne 
are usually expected. The OVII line detected at 0.57 keV 
in the spectrum of MKN 304 is therefore consistent with 
reprocessing in the warm gas. 

However, in sources showing both absorption and emis- 
sion lines the former usually appear to have different veloci- 
ties respect to the latter (and in many cases consistent with 
the galaxy rest-frame, e.g. CoUinge et al. 2001), so it is not 
clear if emitter and absorber systems are intrinsically the 
same. For instance, Kaastra et al. (2002) found the optical 
Narrow Line Region as the most likely origin of OVII and 
Ne IX emission lines observed in NGC 5548. Unfortunately, 
EPIC data do not have enough resolution to shed light on 
this issue. 

We used CLOUDY (Ferland 2001) to further test the 
idea of a possible origin of the oxygen emission line from 
the same warm medium responsible of the absorption in 
MKN 304. We employed the table "AGN" ionizing contin- 
uum, with an X-ray photon index fixed to 1.88, a hydro- 
gen column density set to 9 x 10^^ cm~^, and an ioniza- 
tion parameter equivalent to the highest one obtained from 
XSPEC. For this high ionization state of the gas, the trans- 
mitted X-ray spectrum is only marginally dependent on the 
gas density, the EUV continuum and the UV/X-ray flux ra- 
tio. Following the results of the XSPEC flts, the total output 
spectrum from CLOUDY has been computed as the sum of 
the incident and transmitted spectra with relative contribu- 
tions of 5 and 95%, respectively. In this composite spectrum 
the most prominent emission feature is due to OVII, whose 
EW (~ 110 eV) is in excellent agreement with that inferred 
by our flt, i.e. EW — 105±18 eV. Other emission features in 
the model (i.e. NelX, NVI, MgXI) are several times weaker, 
except for the Fe Ka line, whose EW is similar to that of 
OVII (see Sect. 4.3.). Therefore, we conclude that although 
this result does not ultimately proof the origin of the He-like 
oxygen line in the warm absorber, it clearly supports such 
a scenario. 

Finally, the combination of an unreddened continuum 
in the optical band (Ha/H/3 ~ 3; e.g. Osterbrock 1977) 
and a column density of Nh ~ lO^'^ cm~^ suggests a dust- 
free X-ray absorbing medium in MKN 304, maybe located 
within the dust sublimation radius (see Crenshaw, Kramer 
& George 2003 and references therein). 

4.2 Soft excess 

The soft excess component (whose normalization is « 4% of 
the primary emission) detected in MKN 304 appears similar 
to that observed in many X-ray obscured objects (Turner et 

Nonetheless, it is also worth noting that the colder ionized com- 
ponent has values of Nh , ^ and T approximately consistent with 
the production of UTA features around 0.75 keV. Unfortunately, 
the lack of an adequate energy resolution does not allow us to 
further test this hypothesis. 
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al. 1997) . There are two possibilities to explain the origin of 
such excess, i.e: scattered/reflected emission from the pho- 
toionized gas in the warm absorber and "partial-covering" . 
We discuss in turn these hypothesis in the following. 

XMM-Newton and Chandra results on Seyfert 2 galax- 
ies (Sako et al. 2002 and reference therein) suggest that the 
soft X-ray emission in these objects is the result of primary 
continuum scattering underlying a blend of strong emission 
lines, mostly from H- and He-ions of C, N, O, Ne, and 
Fe. Since warm absorber regions show physical properties 
(i.e. 5, column densities, temperatures) similar to those ob- 
served in the X-ray photoionized emitting plasma in Seyfert 
2 galaxies, it is likely that they are different manifestations 
of the same phenomenon (Kinkhabwala et al. 2002) sim- 
ply seen at different inclination angles. Within this scenario 
the soft excess found in MKN 304 would be therefore due 
to scattered/reflected emission from the ionized outflowing 
plasma. 

Interestingly, Ogle et al. (2000), Bianchi et al. (2003) 
and Iwasawa et al. (2003) detected indeed an extended 
soft X-ray emission dominated by narrow emission lines 
overlapping the optical Narrow Line Region in NGC 4151, 
NGC 5506 and NGC 4388, respectively. According to this 
hypothesis, a significant contribution from this extended soft 
X-ray component to the OVII emission line would be also 
expected (see Sect. 4.1). 

An alternative explanation of this "soft excess" could 
also be some fraction of the primary continuum leaking 
through the clumpy obscuring material (i.e. along the line 
of sight the absorber does not completely cover the nuclear 
source): the so-called "partial-covering" scenario. Further- 
more, as mentioned in Sect. 14.11 dumpiness in the ionized 
absorbing gas appears to be very likely. 

Unfortunately, the present data do not allow to discrim- 
inate these hypothesis. The unabsorbed power law in model 
E (with a slope identical to the primary continuum one) can 
be therefore considered as a good approximation of the real 
(more complex) spectrum. However, given the energy res- 
olution of EPIC CCD cameras and the relatively faintness 
of MKN 304, a description of this soft X-ray component in 
terms of many discrete emission features, as observed in the 
high resolution spectra of some bright AGNs (see Kinkhab- 
wala et al. 2002), cannot be ruled out. 

Finally we rule out the hypothesis of an origin of the soft 
excess from reprocessing (i.e. reflection) by an ionized disk 
since this interpretation is inconsistent with the observed 
centroid of the iron line. In fact, assuming that the line arises 
from such an ionized disk (as expected if a strong refection 
component is present; Ross, Fabian & Young 1999), it should 
be emitted at energy higher than 6.4 keV from iron at high 
ionization states. On the contrary, we found an energy of the 
Fe Ka emission line strictly consistent only with "neutral 
iron" (see below). 

4.3 Iron line 

We report the first clear evidence of an iron fluorescence 
emission line around 6.4 keV in the spectrum of MKN 304. 
The detection of a Fe Ka hue at E = 6.34±0.06 keV con- 
strains the ionization parameter of the emitting medium to 
be 5 5: 100 erg cm~^ s~^, higher values of ^ would indeed 
produce a feature peaked at higher energies (Kallman & 



McCray 1982) excluded by the present data. Interestingly, 
both cold and hot warm absorber components that resulted 
in our best fit (model E) to the EPIC spectrum have values 
of ^ which match well with this constraint. This suggests 
that a significant component of the iron emission line may 
be originated in the warm medium. 

We have run the same CLOUDY model as in Sect. 4.1 
to predict the strength of the Fe Ka emission line. The ex- 
pected value (~ 110 eV) falls in the range we measured from 
the XMM-Newton spectrum (i.e. EW = 79±35 eV). So, in- 
terestingly, the iron line detected in MKN 304 is consistent 
with an origin in the warm absorber observed in this source. 

We have also tried to fit the data with a model in- 
cluding a neutral reflection Compton component (pexrav 
model in XSPEC) in order to test a possible origin of the 
iron line from the optically-thick accretion disk. This spec- 
tral parametrization yielded a good fit with a x^i'^) ~ 
335(357): however, there is no statistical improvement re- 
spect to Model E. Although the detection of cold reflec- 
tion results therefore only marginal, we cannot rule out its 
presence. The observed spectral range (0.3-12 keV) of this 
EPIC observation is not adequate to constrain the strength 
of this "bump-like" reprocessed emission feature peaked at 
« 30-40 keV (Lightman & White 1988). The resulting up- 
per limit for the covering factor of the material irradiated by 
the X-ray source is K =^/2n < 0.7. Both the value of the 
equivalent width of the line and R suggest that the reflec- 
tor - if existing - subtends a solid angle of tt steradian 
to the X-ray source, i.e. intercepts ^ 50% of the primary 
continuum once a slab geometry for the reflecting medium 
is assumed (George & Fabian 1991). 

Furthermore the narrow proflle of the Fe Ka line 
(cTKc <0.18 keV corresponding to Aw(FWHM) < 8000 
km/s for the line velocity width) implies that even if the 
reprocessing medium is the accretion disk, the fluorescence 
does not take place from its innermost (relativistic) regions. 
A substantial contribution to the Fe Ka emission could be 
provided from the putative molecular torus located at 1-10 
pc from the central X-ray source as invoked in the Unifica- 
tion model of AGNs (Antonucci 1993). In this case the ma- 
terial is assumed to be Compton-thick (i.e. with Nh ^ 1.5 
X 10^" cm"^, see Matt, Guainazzi & Maiolino 2003), and 
the emission is due to the refiection of the primary contin- 
uum off the inner walls of the torus into our line-of-sight 
(Ghisellini et al. 1994). 

The observed properties of the Fe Ka line in MKN 304 
appear to be very similar to the results obtained by the 
spectral analysis of the brightest Seyfert 1 galaxies using 
XMM-Newton (Pounds & Reeves 2002) and Chandra (Pad- 
manabhan & Yaqoob 2002) data. These works indeed shows 
the "ubiquitous" presence of a narrow iron line component 
centered at 6.4 keV with an EW « 50-100 eV and axa <120 
eV in the X-ray spectra of the low-luminosity (i.e. with 
1/2-10 < 10*^ erg s-i) AGNs. 



5 SUMMARY 

The ~30 ks XMM-Newton exposure presented here is 
the first X-ray broad-band (0.3-12 keV) observation of 
MKN 304. Our data analysis has revealed that this Seyfert 
1 galaxy has a complex spectrum, dominated by heavy ion- 
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ized obscuration. This warm absorber can be well fit by a 
"two-zone" gas with ionization parameter of ^ ~ 6 erg cm~^ 
s~^ and ^ ~ 90 erg cm~^ s~^, respectively. We estimate that 
the warm absorber has a column density of Nh ~ 9 x 10^^ 
cm~^ in the high-ionization component and Nh ~ 2 x 10^^ 
cm~^ in the low-ionization component. Such values are cur- 
rently amongst the largest ones observed so far in Seyfert 1 
galaxies by XMM-Newton and Chandra. 

A particularly interesting result of our analysis is the 
detection of an emission feature at ~ 0.57 keV, likely due 
to the helium-like oxygen from the warm absorber itself or 
regions much farther out from the central BH (i.e. the optical 
Narrow Line Region). 

We also report the first detection (significant at 99.9%) 
of an iron line in this source. The combination of the energy 
centroid (« 6.4 keV) and the upper limit of the line width 
(cT/fd < 0.18 keV) suggests an origin from "cold" (i.e. Fel- 
XV) matter distant from the inner disk region. Alternatively, 
the line strength is consistent with the emission arising from 
the warm absorber itself. 

Finally, we have also found a weak (4% of the primary 
continuum) soft excess. The presence of this excess could be 
explained by either emission/scattering from a warm gas or 
partial covering, or a combination of them. 
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APPENDIX A: HIGH BACKGROUND 
FILTERING PROCEDURE 

XMM-Newton experiences high "flaring" background pe- 
riods due to the crossing of clouds of high-energy parti- 
cles (mainly soft protons), which are focussed through the 
mirrors to the detectors (Lumb et al. 2002). Time inter- 
vals affected by background flares are essentially useless for 
the scientific analysis. The standard method to fllter out 
these high background periods consists in rejecting peri- 
ods of high count rate (CR) at energies ^ 10 keV (see 
XMM-Newton Science Analysis System Users Guide v2.1; 
Loiseau 2003). The recommended thresholds for the back- 
ground count rate are 0.35 count/s and 1 count/s for the 
MOS and the PN camera, respectively. This approach (i.e. 
using flxed threshold for the CR of the background) well 
adapts for the detection of faint sources but it could produce 
very conservative results in the case of high-flux sources. 

In the reduction of the MKN 304 data, we thus apply 
the following alternative procedure ("MaxSNR" hereafter) 
to reject the high background level periods. The underlying 
idea basically consists of flltering out just those time inter- 
vals for which the count rate of the background reach values 
so high that the SNR of the source does not significantly 
improve (or even worses) by including such time intervals 
in the analysis. To this aim, lightcurves of the source and 
background region, i.e. a proper annulus around the source 
large enough to provide a good estimation of the background 
(with a 120 arcsec radius in the case of MKN 304), are ex- 
tracted in the energy band for which the SNR should to 
be maximized. In the analysis of MKN 304 data, we have 





Standard method 


MaxSNR method 


Final exposure 


6.7 ks 


29 ks 


Obs. countrate 


0.50 ±0.01 


0.480 ±0.005 


SNR 


53 


103 


Total counts 


2971 


16077 


(0.3-12 keV) 







Table Al. Comparison between the results obtained applying 
the "standard" method and the "MaxSNR" method presented in 
this paper for the high background filtering in MKN 304. Data 
are reported for the PN camera. 




Figure Al. Comparison between the PAf spectrum of MKN 304 
obtained using the MaxSNR filtering method (top) and the stan- 
dard filtering procedure {bottom). After the filtering the exposure 
times result of 29 ks and 6.7 ks, respectively. 



taken into account the 0.5-1.3 keV (0.6-1.3 keV) band for 
PN{MOS) in order to optimize the SNR in the band where 
most of the prominent warm absorber features are located. 

On the contrary, the "standard" filtering method re- 
quires to extract a lightcurve of the entire field of view for 
energies above iJ > 10 keV. 

Lightcurve datapoints for the source and background 
are sorted for increasing values of the CR of the background 
region. The cumulative signal-to-noise (SNRcum) distribu- 
tion function is then calculated. 

The SNRcum as a function of the increasing background 
CR builds a curve which monothonically rises up to a peak. 
After this maximum value, the SNRcum distribution indeed 
fiattens or even decreases. This value is the threshold used to 
reject the high background periods. In fact, if included, these 
periods would not provide any significant improvements to 
the final SNR of the spectrum of the source, but it would 
increase the background contribution to the raw source spec- 
trum. 

This method to filter out high background periods op- 
timizes the duration of the observation useful for a mean- 
ingful scientific analysis on the basis of the ratio between 
the count rates of source and background. It appears par- 
ticularly appealing in the case of bright source for which 
an even "flaring" background can represent just a negligi- 
ble fraction of their CR. It also has the advantage that the 
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filtering is performed in a particular energy band (chosen a 
priori) where the SNR need be maximized for the accurate 
study of peculiar spectral features. 

Table IXn lists the final values of exposure, total counts 
in the 0.3-12 keV band, countrate and SNR for the PN 
camera obtained using the "MaxSNR" method compared 
to those derived from the application of the "standard" 
method. In particular, the method presented here yields a 
number of counts which is about 5.5 times larger than that 
resulted from the "standard" method. Fig. lAll shows the 
comparison between the spectra obtained using the two fil- 
tering procedures: the higher quality of the PA'' spectrum 
provided by the "MaxSNR" method is clearly evident. 

This paper has been typeset from a TjjX/ I^T^^X file prepared 
by the author. 



